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Since Bridgman’s discovery of five solid water (HO) ice phases’ 
in 1912, studies on the extraordinary polymorphism of H20 have 
documented more than seventeen crystalline and several amorphous 
ice structures”, as well as rich metastability and kinetic effects*>. 
This unique behaviour is due in part to the geometrical frustration 
of the weak intermolecular hydrogen bonds and the sizeable 
quantum motion of the light hydrogen ions (protons). Particularly 
intriguing is the prediction that HO becomes superionic® !7—with 
liquid-like protons diffusing through the solid lattice of oxygen— 
when subjected to extreme pressures exceeding 100 gigapascals 
and high temperatures above 2,000 kelvin. Numerical simulations 
suggest that the characteristic diffusion of the protons through the 
empty sites of the oxygen solid lattice (1) gives rise to a surprisingly 
high ionic conductivity above 100 Siemens per centimetre, that 
is, almost as high as typical metallic (electronic) conductivity, 
(2) greatly increases the ice melting temperature’~!° to several thousand 
kelvin, and (3) favours new ice structures with a close-packed 
oxygen lattice", Because confining such hot and dense H20 in the 
laboratory is extremely challenging, experimental data are scarce. 
Recent optical measurements along the Hugoniot curve (locus 
of shock states) of water ice VII showed evidence of superionic 
conduction and thermodynamic signatures for melting’, but 
did not confirm the microscopic structure of superionic ice. Here 
we use laser-driven shockwaves to simultaneously compress and 
heat liquid water samples to 100-400 gigapascals and 2,000- 
3,000 kelvin. In situ X-ray diffraction measurements show that 
under these conditions, water solidifies within a few nanoseconds 
into nanometre-sized ice grains that exhibit unambiguous evidence 
for the crystalline oxygen lattice of superionic water ice. The X-ray 
diffraction data also allow us to document the compressibility of 
ice at these extreme conditions and a temperature- and pressure- 
induced phase transformation from a body-centred-cubic ice phase 
(probably ice X) to a novel face-centred-cubic, superionic ice phase, 
which we name ice XVIII”. 

Our experimental system (Fig. 1) builds on the development of 
advanced dynamic compression experiments to access high-density, 
high-temperature states beyond the reach of current static compres- 
sion methods, and uses the nanosecond, in situ, powder X-ray diffrac- 
tion (XRD) techniques at the Omega Laser facility'®'°. Water being 
quite compressible, single-shock temperature rises very quickly with 
pressure’® reaching 6,000 K near 100 GPa (Fig. 1d). Therefore, to cre- 
ate ultrahigh-pressure solid ices in the laboratory, a series of multiple 
shocks is necessary to realize a quasi-isentropic compression and to 
minimize the temperature increase”™”!. 

Six high-power laser beams are fired in a temporally tailored 
sequence about 15 ns long to launch a series of compression shock- 
waves into a thin water layer encapsulated between a diamond ablator 
and a diamond window (Fig. 1a). The multiple shockwaves reverberate 
between the two stiff diamonds to achieve a homogeneous quasi- 
isentropic compression of the water layer?” (Fig. 2a). 


Adjusting the energy and timing of the drive-laser beams on each 
laser shot, water is compressed to ten different pressure-temperature 
states in the ranges 100-420 GPa and about 2,000-3,000 K. Ultrafast 
Doppler velocimetry monitors the compression history on each shot: 
a two-channel, line-imaging Velocity Interferometer System for Any 
Reflector (VISAR) records the arrival and exit times of the leading 
shock in the water sample layer, and encodes the time history of the 
velocity of the diamond window free surface, ups(t), which is progres- 
sively accelerated by the arrival of the successive compression shock- 
waves (Figs. 1c, 2a, and Supplementary Fig. 1). 

Radiation-hydrodynamic simulations are performed for each 
laser experiment (Fig. 2). An iterative optimization of the drive-laser 
intensity input allows us to reproduce the measured timing and ampli- 
tude of the successive shocks determined by analysis of the velocime- 
try data. Once a good match between the measured and computed 
free-surface velocity history ups(t) is obtained (Supplementary Fig. 3a), 
we interrogate the simulations to determine the pressure-temperature 
compression path experienced by the H20 sample and the pressure- 
temperature conditions during the XRD snapshot using the timing of 
the X-ray pulse (Supplementary Fig. 3). The simulations indicate that 
the XRD data were collected along two quasi-isentropic paths near 
2,000 K and 3,000 K, depending on the strength of the first shockwave 
reaching the water layer and on the timing of the main compression 
pulse. 

The atomic structure of dynamically compressed H20 is character- 
ized in situ using angular dispersive XRD recorded on multiple image 
plate detector panels (Fig. 1b). Focusing 16 laser beams on an iron 
(Fe) foil generates an intense X-ray pulse at 6.7 keV resulting from 
helium (He)-like atomic transitions. About 2 x 10!” X-ray photons 
are collimated by the 500-j1m-diameter tantalum (Ta) aperture placed 
immediately behind the target package (Fig. 1a). Diffraction off the 
walls of the Ta pinhole imprints a zero-pressure XRD reference pattern, 
superimposed on the XRD signal coming from the H2O sample on the 
image plates. On each shot, the timing of the 1-ns-long X-ray pulse is 
adjusted to capture the structure of water near peak compression, typ- 
ically 10-15 ns after the first shockwave arrival in H20 (Figs. 1 and 2a). 

Figure 3 shows representative XRD images from four different exper- 
iments, projected into the Bragg angle (20)—azimuthal angle (4) space. 
One diffraction line, which we attribute to H20 ice, is found in the XRD 
patterns recorded for all ten shots and demonstrates the solidification of 
a macroscopic fraction of the water sample into a crystalline solid (see 
Supplementary Information sections S4 and S5 for a detailed descrip- 
tion of all the XRD data). 

Interrogating the hydrodynamic simulations reveals that the growth 
time delay (At) between the instant H20 is brought into the stability 
field of ice and the time of the XRD snapshot (see Fig. 2b and values 
in Supplementary Table 2) ranges between 2 ns and 5 ns in our experi- 
ments. This puts an upper bound on the nucleation of ice near 2 ns. This 
is faster than previous studies on water ice at much lower pressure?” 
and can be interpreted as a signature for the large over-pressurization 
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Fig. 1 | Experimental concept (not to scale). a, A series of laser pulses 
generates a multi-shock compression sequence of the initially liquid 
water layer. Additional laser beams generate a X-ray flash off an iron foil. 
b, Image plates record the X-ray diffraction snapshot documenting the 
solidification of the initially liquid water during the dynamic loading. 
Ambient-pressure diffraction lines from the Ta pinhole are used to 


(far into the solid stability field): the large free-energy difference favour- 
ing ice over liquid water drives the rapid solidification. 

Figures 1b and 3a also reveal that the diffraction peaks from H2O 
are broader than the zero-pressure Ta reference lines, whose width is 
dominated by the experimental angular resolution in 20 and geomet- 
rical effects (contributing to about 0.8°). Because we observe only the 
most intense XRD peak for water ice, we cannot formally disentangle 
the respective strain and grain-size contributions to the peak broad- 
ening. However, assuming that grain-size effects solely determine the 
angular broadening and using the Scherrer equation with the full- 
width at half-maximum (FWHM) of the H20 line provides a lower 
bound for the typical grain (crystallite) size D of about 4-7 nm (see 
Supplementary Information sections S4 and S6). This is similar to 
recent observations in laser-driven shock crystallization of glassy SiO3 
(ref. *4) and is also consistent with billion-atom molecular dynamics 
simulations of strongly undercooled iron melt, indicating freezing by 
homogeneous nucleation within about 200 ps followed by competitive 
coalescence growth” of approximately 5-10-nm-sized grains within 
1-2 ns. 

The d-spacing of the observed H2O XRD line and its decrease upon 
pressure increase from 101 GPa to 249 GPa along the quasi-isentrope 
near 2,000 K are consistent with the most intense (110) reflection of the 
body-centred-cubic (b.c.c.) oxygen lattice of ice X (Fig. 3b, blue circles). 
Indeed, the experimental d-spacing and density—pressure dependences 
agree with both experimental-based”® and density functional theory 
(DFT)-based equation-of-state models?’ (black dashed line and blue 
dotted line in Fig. 3b and c, respectively). 

At higher pressure along the same quasi-isentrope near 2,000 K, 
the data reveal a jump in the d-spacing of about 0.05 A between 
249 GPa and 278 GPa followed by a compression-driven decrease up 
to 420 GPa, pointing to a structural phase transformation (red circles 
in Fig. 3b). The observed XRD line above 278 GPa most probably cor- 
responds to the (111) reflection off a close-packed ice structure with a 
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calibrate the diffraction geometry (data from shot 73144). c, A line- 
imaging velocimeter (VISAR) monitors the diamond window free-surface 
velocity history ups(t) from which the pressure-temperature compression 
path is reconstructed. d, The reverberation compression scheme enables 
the phase diagram of H20 at extreme pressure-temperature conditions to 
be explored. 


face-centred-cubic (f.c.c.) lattice of oxygen atoms (see Supplementary 
Table 2). Attributing the observed XRD line to other potentially 
expected structures would result in very large, unrealistic density 
jumps above 249 GPa (see Supplementary Information section S5). 
Interpreting the new phase as f.c.c. does not produce a substantial dis- 
continuity in the pressure—density relation for dense H20 at the b.c.c.- 
f.c.c. transition (within the uncertainty of our measurement), consistent 
with the expected small volume jump for this phase change", indicative 
of the novel f.c.c. structure having a volume compressibility similar to 
that of ice X. 

The two XRD patterns recorded near 3,000 K (red squares in Fig. 3b, 
c) are also consistent with the same f.c.c. structure. Showing the results 
of the XRD experiments on the phase diagram of dense HO (Fig. 4) 
indicates that both pressure (between 249 GPa and 278 GPa) and tem- 
perature (between about 2,000 K and about 3,000 K) can induce the 
transformation of b.c.c. H2O ice into the novel f.c.c. phase that we pro- 
pose to name ice XVIII (see Supplementary Table 1 in Supplementary 
Information section S7). 

Figure 4 also shows that the pressure and temperature of (at least) 
one of our XRD measurements—near 160 GPa and 3,200 K—overlap 
with the conditions where superionic conduction has recently been 
experimentally observed!6. The new ice (ice XVIII) is therefore very 
likely to be superionic and the XRD data provide direct experimental 
evidence for a crystalline sublattice formed by the oxygen atoms in 
superionic water ice. 

The present XRD dataset also confirms two additional predictions 
from computer simulations for superionic water ice. First, it confirms 
the prediction’?! that the extra entropy associated with the mobile 
protons favours novel close-packed structures—such as f.c.c. ice 
XVIII—above about 100 GPa. In contrast, a f.c.c. lattice is not expected 
in insulating ice in this pressure range because body-centred-cubic ice X 
should be stable to at least 300 GPa, before transforming to space group 
Pbcm?8-3!, Second, it confirms the prediction®"'>!>!> that superionic 
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Fig. 2 | Determination of the reverberation compression path using 
one-dimensional Lagrangian radiation hydrodynamic simulations. 
a, Example space-time contour map of the longitudinal stress for shot 
75758: a series of shockwaves compresses the initially liquid water layer 
before the XRD snapshot at t ~ 11 ns (white rectangle). b, Calculated 
pressure-temperature compression paths experienced by each simulation 
mesh element (dotted lines), and average for the water layer (solid 
line), with time as a colour scale. The XRD snapshot was taken just 

a few nanoseconds after H2O has been dynamically compressed into 
the predicted stability field of superionic ice!" (black arrows). The 
final estimate of pressure and temperature during the XRD snapshot is 
indicated by the blue square, with error bars representing the standard 
deviation of the distribution of states during the X-ray exposure. 


ice can sustain higher temperatures before melting at high pressures, 
leading to a rapid increase of the melting line above 1,000 K and 50 GPa. 
Our XRD data up to 320 GPa and 3,800 K provide unambiguous lower 
bounds for the melting line of water ice at high pressure (Fig. 4), in 
agreement with previous static compression experiments at lower pres- 
sure and the observation of two thermodynamic anomalies near 4,800 K 
and 190 GPa along the Hugoniot of pre-compressed ice VII’®. 

Figure 4 further indicates that the pressure-temperature conditions 
for the b.c.c. phase data lie in the vicinity of the predicted boundary 
between solid insulating and superionic ices'® (grey dashed line). Since 
temperature-driven phase transitions are more likely to favour trans- 
formations from a close-packed to an open, higher-entropy structure 
upon temperature increase (for example, f.c.c. to b.c.c.), rather than 
b.c.c. to fic.c. (as we observe in our experiments), the data collected in 
the b.c.c. phase are better interpreted as insulating solid ice. If, from 
a thermodynamic point of view, a temperature-induced superionic 
b.c.c.-to-superionic f.c.c. transition seems unlikely, the extra entropy 
of the diffusing protons could, on the other hand, favour a tempera- 
ture-driven transformation of insulating b.c.c. to superionic f.c.c. near 
200 GPa. 
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Fig. 3 | X-ray diffraction data and results. a, Representative data projected 
onto 20-4 space and the overlaid one-dimensional integration showing the 
evolution of the diffraction line for water (asterisks) with pressure along 
the quasi-isentrope near 2,000 K (blue for b.c.c. and red for f.c.c.). See 
Supplementary Information section S4 and Supplementary Figs 5 and 6 for 
the complete peak assignment. b, Pressure dependence of the measured 
d-spacing. Literature data are taken from refs 7°’. A discontinuity is 
observed at 278 GPa along the quasi-isentrope near 2,000 K (blue and 

red circles for b.c.c. and f.c.c. respectively). Near 3,000 K the d-spacing is 
consistent with the high-pressure f.c.c. structure (red squares). Typical 
uncertainties are less than about 0.01 A, smaller than the size of the 
markers. DAC, diamond anvil cell. c, Pressure-density experimental data 
compared with previous data and computed isotherms. Data are taken from 
refs 1526273337, Error bars represent lo uncertainty for each experiment; 
see Supplementary Information for a detailed discussion. 


The XRD data collected near 100 GPa and 2,000 K—that is, along the 
ice VII Hugoniot, very close to the conditions where superionic con- 
duction was experimentally documented '°—might indicate that a b.c.c. 
superionic ice exists below 160 GPa, as suggested by simulations” =". 
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Fig. 4 | Phase diagram of H20. Pressure-temperature conditions of the 
XRD measurements in the b.c.c. phase (blue circles) and in the novel 
f.c.c. superionic ice, ice XVIII (red circles and squares, near 2,000 K and 
3,000 K, respectively). Conductivity measurements with precompressed 
shocks'® (magenta squares) and reverberating shocks”! (magenta 
diamonds) have demonstrated superionic conduction between 100 GPa 
and 190 GPa along the ice VII Hugoniot. Coloured crosses indicate 


However, more data at higher temperature in this pressure range are 
needed to remove any ambiguities. 

This study brings the last missing piece to the puzzle regarding the 
existence of superionic water ice, 30 years after the original predic- 
tion®, and directly documents the presence of a crystalline lattice in 
superionic ice at the extreme pressures and temperatures relevant 
for the deep interior of Neptune, Uranus and icy giant exoplanets. 
Given that H20 is believed to constitute about 65% by mass of Uranus 
and Neptune, comparing the interior isentropes from ref. ** with the 
experimentally constrained melting line on Fig. 4 confirms the predic- 
tion!™! b13 that superionic ice could indeed be stable at depths greater 
than about 8,000 km (roughly one-third of the radius of Uranus and 
Neptune) and contribute to forming a thick solid icy mantle in these 
planets. 

The observation of a crystalline lattice of oxygens indicates that 
superionic water ice should indeed behave similarly to a solid, with a 
viscosity several orders of magnitude larger than typical fluid values. 
This finding is supported by computing the elastic constants with DFT- 
based molecular dynamics (DFT-MD) up to 450 GPa at 4,000 K (see 
Supplementary Information section S8). The calculation shows that the 
shear modulus of superionic ice is slightly lower than that of insulating 
ice X but evolves similarly with increasing pressure (see Supplementary 
Fig. 24). Therefore, the characteristics of superionic ice flows within 
Neptune and Uranus are likely to be more similar to those of slow con- 
vection within the solid Earth mantle than to those of the rapid swirls 
of fluid iron in the Earth’s outer core. 
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pressure-temperature conditions for static compression XRD experiments 
(see details in Supplementary Fig. 18). Planetary interior conditions for 
Neptune and Uranus from ref. ° are also shown. The predicted boundary 
between solid and f.c.c. superionic ice" is represented by the dashed 

grey line. Error bars represent 1o uncertainty for each experiment; see 
Supplementary Information for a detailed discussion. 
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